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Abstract 



For the tendering of long-term transportation contracts in the bulk industry, ship- 
pers use bidbooks that specify for each lane the load location, the destination, the 
product and the volume that has to be transported over the next so many years. 
Bidbooks are sent out to a preselected group of carriers, who subsequently quote a 
price for each of the lanes. After the return of the bidbooks, the shipper determines 
the winning carriers. Although price is the main driver, a shipper may take into 
account many additional considerations, including upper bounds on the number 
of winning carriers in total, per load location, per country of destination, and the 
maximal transport volume any carrier is allowed to win. The winner determination 
problem is the problem of finding an allocation of the lanes to the carriers so as 
to minimize total transportation casts such that each lane is assigned to exactly 
one carrier and the additional constraints are satisfied. The winner determination 
problem is extremely hard in practice, and most carriers resort to simple rules of 
thumb to select the winning carriers. This is no surprise: we prove that the problem 
is Jv7>-hard in the strong sense. We model the winner determination problem as an 
integer linear programming (ILP) problem and try and solve the model by use of 
CPLEX, a state-of-the-art ILP solver. Tested on a comprehensive set of instances 
generated along the characteristics of a rea.l world case of a European chemical 
shipper with about 4,000 lanes, CPLEX is capable of solving instances with no 
more than 270 lanes to optimally, underlining the practical difficulty of the winner 
determination problem. However, we also develop a fast randomized heuristic, and 
we show that it performs remarkably well, with a gap of no more than 0.8% from 
optimality. This performance supports our conclusion that our heuristic can be 
used to obtain good approximations quickly for realistic, much larger instances. 



1 Introduction 

For tendering long-term transportation contracts in the bulk industry, shippers use 
bidbooks that specify for each lane the load location, destination, the product and 
the volume that has to be transported. The bidbooks are sent out to a preselected 
group of carriers that meet a certain set of minimal quality requirements. These 
carriers then quote a price for each of the lanes. After the return of the bidbooks, 
the shipper determines the winning carriers so as to minimize total transportation 
costs. 

Although price is the main driver in this decision, it is not the only aspect that 
plays a role. A shipper typically takes into account many additional considera- 
tions, as an otherwise inefficient or undesirable situation from an operational and 
transaction cost point of view may occur. These include an upper bound on the 
number of winning carriers in total, per load location, per country of destination, 
and the maximal transport volume any carrier is allowed to win. The problem of 
allocating lanes to carriers in order to minimize transportation costs subject to such 
considerations is called the winner detcmination problem. 

Price is a major driver in the tendering, because transportation of chemical 
bulk products is basically a commodity, or in Kraljic (1983)'s terminology, a leverage 
product, which means that the financial impact on the shipping company is high and 
the risk of not finding a carrier is low. To illustrate this, the last tri-annual tendering 
of the European chemical company that motivated our research concerned about 
€100,000,000, whereas the size of its group of preselected carriers was about 30. 
The role of price makes e-marketplaces and reverse auctions particularly suitable 
for the sourcing of such transportation services; see for instance Van Weele (2000) 
for the impact of Internet technology on procurement. 

Currently, many transportation c-marketplaces exist, in which shippers and 
carriers participate to match orders for and bids on transportation services. Ship- 
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pers participate to find a carrier for certain transportation services and carriers 
participate to find freight for an otherwise empty return trip. However, most exist- 
ing transportation marketplaces are used for spot sourciny, that is for short term 
transportation contracts, see Kaplan and Sawhney (2000). Examples of indepen- 
dent transportation e-marketplaces are CargofirKler.com, Logistics.com, and Open- 
ship.com. Most of these marketplaces do not offer services to facilitate shippers' 
systematic sourcing, that is long term transportation contracts. A notable excep- 
tion is Translogistica.com. In our paper, we focus on the systematic sourcing of 
transportation services. 

The procurement of transportation contracts is extensively described by Caplice 
(1996), who also explores the possibilities for operating a reverse auction for ten- 
dering transportation services. Furthermore, he presents various carrier assignment 
models that include some of the considerations we pointed out earlier. Ledyard et al. 
(2002) describe a so-called combined- value auction, that is, a combinatorial auction 
for transportation services for Sears, Roebuck and Co., hut hardly no attention is 
paid to modeling and algorithmic aspects. In business practice, Logistics.com and 
CombineNet.com have shown that an auction of transportation services of this type 
is really beneficial. 

Our work is motivated by a European chemical company that explored the 
possibilities of having an electronic reverse auction for tendering transportation 
services. Just a virtual bidbook on the Internet where carriers can submit their 
bids ; is not sufficient to operate such a reverse auction. Other issues that have to 
be addressed before operating such an auction include the design of the auction 
mechanism (for instance, single round versus multiple-round; see Cramton (1998) 
and the references therein) and the solution of the winner determination problem. 

Approximating rather than optimizing the winner determination problem has 
consequences for the economic efficiency of the auction. Nisan and Ronen (2000) 
show that under a VCG-inechanism a combinatorial auction is truthful if and only 
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if the winner determination problem is solved to optimally. A truthful mechanism 
maximizes the total welfare, and hence approximating the winner determination 
problem can lead to a lower social welfare and revenue loss. 

Accordingly, there is a growing body of literature on optimization algorithms 
for the winner determination problem, in particular for so-called combinatorial auc- 
tions, where bidders are allowed to bid XOR bids (I want X or Y but not both). 
Sandholm (2002) developed a search algorithm to solve the winner determination 
in a combinatorial auction to optimally. The search generates each allocation 
with positive revenue exactly once, which results in a worst case complexity of 
C(n m ), where n is the number of bids and m the number of items on sale. Fu- 
jishima, Leyton-Brown, and Shoham (1999) present a depth-first search algorithm 
with a number of speed-ups for the winner determination problem in combinato- 
rial auctions. Davenport and KaJagnanam (2002) present a volume discount and a 
combinatorial auction mechanism for the procurement of direct inputs of a manu- 
facturing firm. For both auctions they present a mathematical programming model 
of the winner determination problem. The models take into account several of the 
side constraints that are also present in our problem, such as the maximum number 
of suppliers and the maximum quantity to be allocated to any supplier. Andersson 
et al. (2000) present a mixed integer linear programming formulation for this win- 
ner determination problem, which can be solved by standard mixed integer linear 
programming solvers. They compare this approach with existing algorithms on dif- 
ferent problem distributions and conclude that, such a standard approach performs 
quite reasonably. Bikhchandani and Ostroy (2001) present a number of different 
extended linear programming formulations for the winner determination problem, 
one of which always has an integral solution. Bikhchandani et ah (2002) propose an 
other extended integral formulation with fewer variables. Nisan (2000) suggests a 
linear programming-based approach for solving the winner determination problem 
under different bidding languages. The method first solves a linear program, af- 
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Ler which branch-and-bound is applied to find a feasible assignment. Giinluk et al. 
(2002) present an alternative integer linear programming formulation based on a set 
packing formulation with a larger number of variables, which gives rise to a tighter 
linear programming relaxation. In this formulation, variables indicate whether a 
bundle is assigned to a bidder or not. The authors present also a branch-and-price 
algorithm, which served as shadow-solver in one of the FCC-auctions. Hoos and 
Boutilicr (2000) present a stochastic local search algorithm for the winner determi- 
nation problem and show that this method finds high quality solution much faster 
than current optimal methods. Zurel and Nisan (2001) present an approximation 
algorithm for the winner determination problem in a combinatorial auction that 
consists of approximating the linear programming relaxation and then applying a 
sequence of greedy hill-climbing algorithms to improve on the initial solution. The 
authors claim an average approximation error of less than 1% for a variety of prob- 
lem instances. For an extensive overview of the winner determination problem in 
combinatorial auctions, we refer to De Vries and Vohra (2003). 

In our problem, the bids are OR-bids for single lanes, and not XOR. The 
comphcating feature of our problem lies in the many combinatorial constraints on 
the possible allocations to the carriers. These two aspects differentiate our problem 
from the existing literature. Furthermore, we have access to the real data of a very 
large tendering of transportation services. 

Indeed, the winner determination problem is extremely hard in practice, and 
most carriers resort to simple rules of thumb to select the whming carriers. This 
is not surprising: the problem is Afp-hard in the strong sense. We model the 
winner determination problem with its many practical side coastraints as an integer 
linear programming (ILP) problem. Tested on a comprehensive set of instances 
generated along the characteristics of a real world case of a European chemical 
shipper with about 4,000 lanes, CPLEX ; a state-of-the-art ILP solver, is capable 
of solving instances with no more than 270 lanes to optimahty, underlining the 
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practical difficulty of the winner determination problem. However, we also develop 
a fast randomized heuristic, and we show that it performs remarkably well. 

The theoretical and practical relevance of our findings is multifold. First, we 
prove that the winner determination problem with OR-bids of single items together 
with many side-constraints is AfP-hard in the strong sense and show also that gen- 
era] optimization technology (such as CPLEX) is unable to .solve the very large 
instances of the winner determination problem to optimally; accordingly, if we 
really wish to apply such technology, then we need to reduce the instance size by 
bundling certain lanes (like Ledyard et ah (2002) have done). Second, the linear 
programming relaxation of our model formulation gives very strong lower bounds, 
- on average no more than 0.6% away from optimal. Third, a reasonably simple ran- 
domized heuristic performs remarkably well, within about 0.8% from optimally. 
Hence, an alternative approach to tackle large, realistic instances would be to use 
the presented heuristic, to run it more often, or to develop a more sophisticated 
version. Finally, we have shown that substantial savings are possible using algorith- 
mic support for the winner determination problem; this confirms similar findings 
by Ledyard et al. (2002) for the Sears, Roebuck and Co. case. 

The outline of the paper is as follows. In Section 2 we describe the problem. In 
Section 3 we formulate the winner determination problem as an integer linear pro- 
gramming problem. Section 4 describes our randomized heuristic. Computational 
results are discussed in Section 5. Section 6 gives the conclusions of the paper. 

2 Problem description 

We consider a transportation tendering process, in which each lane is described by 
the load location, the destination, and the product and volume that have to be 
transported over a pre-specified period of time. Assume T transportation compa- 
nies, referred to as carriers, are bidding on the L lanes. The lanes have desiii.a- 



5 



tions in C different countries and have to be transported from O origins, called 
load locations. The total volume of the products to be transported is v and let 
Vi be the volume that has to be transported for lane i (i = 1, . . . , L). Let 6 0 - 
(i = \,...,L,j = l,...,r) be the bid of carrier j on lane t. We assume that all 
bids are positive integers. If carrier j did not bid on lane t, then 6y = oo. The 
shipper wants no more than MC winning carriers in total. Analogously, no more 
than c fc carriers are allowed to win lanes to country k (k = 1, . . . . C), and no more 
than Ph carriers from load location h> (h = 1, . . . , 0). In order not to become too 
dependent on a single carrier, no carrier may receive a volume of more than the 
fraction q of the total transport volume v. 

We assume that each carrier has sufficient transportation capacity for the lanes 
it wins. If needed, capacity constraints for each carrier could easily be included 
in the model. A feasible allocation is one such that every lane is allocated to 
exactly one carrier, each carrier receives a set of lanes with a total volume no 
more than qv, the number of winning carriers is at most MC, for each country k 
(fc = 1, . . . , C) there are no more than c k winning carriers, and for each load location 
h (h = 1, . . . , O) there arc no more than p h winning carriers. Depending on the 
model parameters it is possible that no feasible allocation exists. If no feasible 
allocation exists, the shipper typically would increase A/C, c k or p h for some k or h 
to enlarge the feasible region. In a typical tendering process, the shipper plays with 
these parameters to evaluate different, scenarios. The objective is to find a feasible 
allocation with minimal transportation costs. 

3 Mathematical formulation 

We formulate the winner determination problem as an integer linear programming 
problem. Let z {j be a binary variable that adopts the value 1 if carrier j is allocated 
lane i and the value 0 otherwise. Let Vj = 1 if carrier j is allocated at least one 
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lane; Vj = 0, otherwise, c* = 1 if carrier j gets at least one lane to country k and 
cf = 0, otherwise. Let = 1 if carrier j gets at least one lane from load location 
n and — 0, otherwise. 

The goal is to find an allocation of lanes to carriers such that every lane is 
allocated to exactly one carrier, all constraints are satisfied and transportation 
costs are minimized. Mathematically, the winner determination problem, in the 
remainder referred to as problem WD, is to find values for the decision variables 
Vj ) £ ij j Pj : and Cj that minimize 



minimize 



r l 

J=ri 1 = 1 



subject to zij-y, <0 *= l,...,L,j = l,...,r (i) 

(2) 



T 



J2vj < MC 

; = 1 
T 

J^Zij > 1 1=1, 
j=l 
L 

Y, v ^i $ W 7 = 1,...,T, 



(3) 



(4) 



OikZij < c) t=l L,j = l,.'.,T,A: = 1,...,C (5) 



T 



(6) 

h in Zij < p 1 ] i = l,...,L,j = l,...,T,h=l,...,0 (7) 

T 

Ep' S P/, A = 1,...,0, (g) 

*y.%ic}\pj l e {0,1} t = i,. = i,...,r. (9) 

Constraints (1) ensure that lanes are only allocated to one of the MC winning 
carriers. Constraint (2) ensures that there are no more than MC winning carriers. 
Constraints (3) guarantee that each lane is allocated at least once. Constraints 



7 



(4) enforce that each carrier is allowed to transport no more than qv in total. 
Constraints (5) ensure that lanes to country k cannot be allocated to carriers that 
are not allowed to transport to country k. Constraints (7) are similar constraints for 
the load locations. Constraints (6) and (8) enforce that there are no more than c k 
winning carriers for each country k (k = 1 , . . . : C), and no more than p h carriers for 
each load location h (h = 1, . . . : Q). Constraints (9) are the integrality constraints, 
ensuring that lanes are not split over different carriers. 

The linear programming relaxation is obtained by relaxing constraints (9). The 
solution of the linear programming relaxation is a lower bound on the optimal 
solution value of WD. 

Next, we prove that the winner determination problem is Af7>-hard in the strong 



sense. 



Theorem 1. The winner determination problem WD is MV-hard in the strong 



sense. 



Proof. The decision variant DWD of WD is the problem: Does there exist an 
allocation of the lanes to the carriers satisfying constraints (1) to (9)? If DWD is 
^-complete in the strong sense, WD is JW-hard in the strong sense. To prove 
that DWD is AfP-complete in the strong sense, we present a polynomial reduction 
from 3-PARTITION, which is known to be A^-complcte in the strong sense (Garey 
and Johnson (1979)), to DWD. The problem 3-PARTITION is denned as follows 

Instance: A set Q = (q u . . . , 93m ) 0 f 3m elements, a bound BgZ+ and a 
size Hq e Z+ for each qeQ such that f < < f and such that £ fl€Q Sq = 
rnB 

Question: Can Q be partitioned into m disjoint sets <?],..., Q m such that 
^qeQi s 9 ~ B f or each t=l,... ,m? 
It is clear that DWD G NV. We will prove that an instance of DWD is a yes- 
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instance if and only if the corresponding instance of 3-PARTITION is a yes-instance. 

For any given instance of 3-PARTITION we construct the following instance of 
DWD 



T 




m, 




L 




3T, 




V 




niB, 




q 




1 

m 




Vi 




s iz for i = 1, . 








1, fori = 1,.. 








T, lbr*= 1,. 




Ph 




T, for h= 1>. 




MC 




T } 




k 




3m. 





Note that ^ < ^ < f . It is obvious that the transformation is polynomial. 
Tn the winner determination problem instance we have m carriers, 3m lanes, each 
carrier can be allocated a volume of maximal B. Each carrier is allowed to transport 
to each country and from each load location. 

Consider an arbitrary yes-instance of 3-PARTITION. Hence, the set Q of el- 
ements can be partitioned into m disjoint subsets Qi,...,Q m such that for i = 
3 , . . . , m holds J^ qeQt s q - B. For the corresponding instance of DWD this means 
that the 'ST lanes can be divided into T disjoint sets A U ...,A T such that for each 
i = 1 , . . . . T holds £j. €4 . vj = qv. Hence, we can assign to every carrier a bundle of 
three lanes with volume qv. This assignment clearly satisfies all constraints. This 
proves that each yes-instance of 3-PARTITION corresponds to a yes-instance of 
DWD. 
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Next, we prove the other implication of the equivalence. Consider a yes-instance 
of DWD that is obtained from a 3- PARTITION instance. In that case, there is 
an assignment satisfying all constraints. Each carrier gets a volume of maximal 
qv. Because mqv = „, he gets exactly qv. Because v { > f , there cannot be four 
lanes or more in a bundle. Since there are m carriers and 3m lanes, every carrier 
gets exactly three lanes. This assignment can be transformed to a solution of the 
corresponding 3-PARTITION problem. Q 

4 Randomized heuristic 

In this section we develop a simple randomized heuristic to find a good feasible 
solution. Note that without loss of generality, we may assume that every carrier 
bids on each lane. After all, if a certain carrier j did not bid on lane t, then we 
simply let 6 0 - = cc. This assumption makes it easier to find feasible solutions. 

We apply the heuristic a pre-specified number of times. The heuristic is as 
follows 

Step 1 Randomly select a carrier. Assign to this carrier each lane for which its 
bid is lowest among all carriers, if the volume of the bundle after adding 
the lane is smaller than or equal to qv and constraints (6) and (8) are not 
violated by the assignment. 

Step 2 If all lanes have been assigned: STOP. 

Step 3 If less than MC carriers have been selected: go to step 1. 

Step 4 For each unassigned lane, determine the set of selected carriers to which 
this lane can be allocated without violating any of the constraints. If this 
set is empty, then abort the heuristic - a feasible solution cannot be found. 
Otherwise, allocate the lane to the carrier with the lowest bid among the 
ones that can be allocated the lane without violating the constraints. 



10 



In theory, it is AfV-hard to find a feasible allocation as we have shown in 
Theorem 1. However, in practice the heuristic is very fast and if q, c k {k = 1, . C) 
and p h (h = 1, . . . , O) are reasonable large, it is likely that a feasible allocation will 
be found. Note that the heuristic is likely to give a different allocation each time 
it is run, because of its randomized nature. 



5 Computational results 

The algorithm was coded in C++. The integer linear programs were solved by the 
software package CPLEX, (ILOG, 2001), version 7.1. The tests were performed on 
a 1000 Mhz PC with 200 MB RAM. 

Wc tested the empirical performance of the integer linear programming ap- 
proach and the heuristic on a comprehensive set of instances drawn from h real-life 
case for a. large European chemical shipper. In this data set, there arc 10 carriers, 
destinations in 27 countries, 10 load locations, 4,000 lanes with the volume per lane 
between 10 and 5,000 tons, and the price per ton for each lane between €20 and 
€500 depending on the distance between the load location and the destination and 
the type of product that has to be transported. Since the market for transportation 
of chemical bulk products is very competitive, bids for a lane do not differ much. 

We used the following method to generate random instances of the winner 
determination problem 

Step 1 For each lane i the volume v { is randomly generated uniformly from the 
integers between 10 and 5,000. 

Step 2 For each lane i the average price per ton r, is randomly generated uniformly 
from the integers between 20 and 500. 

Step 3 The price per ton for lane i by carrier j is randomly generated uniformly 
from the integers between [^fj and [-^J. The total bid of carrier ; for 
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lane i equals this price multiplied by v,. 

For all instances, the maximal fraction q of the transport volume is O.C and 
the maximal number of carriers for each load location and country is 2 (c k = 2, 
k - 1, . . . , C and p h = 2, h = 1 , . . . , O). Instances are generated for MC = 3, 4 
and 5 and T = 5 and T = 10. Instances with 25 countries and 2 load locations, 
with 20 countries and 3 load locations and with 27 countries and 10 load locations 
are considered. For each set of parameters we generated 25 instances. 

We report the average and maximal relative gap between the best found heuris- 
tic solution and the optimal linear programming solution, the average and maximal 
relative gap between the best found heuristic solution and the optimal integer linear 
programming solution, the number of instances out of 25 for which the optimal lin- 
ear programming relaxation solution was integer, the average and maximal compu- 
tation time, the average and maximal time needed to solve the linear programming 
relaxation and the average and maximal number of nodes in the branch-and-bound 
tree; see Table 1 for a summary of our notation. The computational results are given 
in Table 2. For each set of instances the optimal. solution of the linear programming 
relaxation is found within 30 seconds on average. 

The instances with ten carriers take on average about twice as much computa- 
tion time as the instances with five carriers. In part, the difference in computation 
times can be explained by 28% of the instances with five carriers having integral 
solutions, against only 9% of the instances with ten carriers. We found no direct 
relationship between MC, the number of carriers allowed, and the difficulty of an 
instance. For the instances with C = 27, L = 10 and T = 5, the heuristic gives a 
solution equal to the optimal integer linear programming solution for 18 instances 
if MC = 3, for 19 out of the 25 instances if MC = 4, and for 18 out of the 25 
instances if MC = 5. So, our heuristic succeeds in finding an optimal solution quite 
often, which fuels our belief that our heuristic might even work better for larger 



instances. 
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We compute the relative gap as the difference between the heuristic solution 
value and the linear programming lower bound divided by the linear programming 
lower bound. For our heuristic, this gap was only 0.8% on average with an average 
computation time of less than 2 seconds. 

We also ran CPLEX with the solution found by our heuristic as an initial 
solution for the instances with C = 27, O = 10 and T - 5. These results are given 
in Table 3. The average total running times are somewhat smaller; for the smaller 
instances, the improvement is even less. 

6 Conclusion and future research 

Sodhi (2001) stresses the importance of Operations Research for improving the 
efficiency of supply chains, as it can be used in functionalities for facilitating pro- 
curement and planning decisions which can be offered by busincss-to-business e- 
marketplaces additional to the spot sourcing services. Kt^kinocak and Tayur (2001) 
and Geoffrion and Krishnan (2001) point out that algorithmic support is needed to 
match orders and bids in e-marketplaces. 

Tn this paper, we have tried to contribute to this area by addressing a winner 
determination problem in a reverse, auction (without XOR bids) for the tendering 
of transportation services with many practical and complicating side-constraints. 
The disappointing performance of the optimization method stresses that still many 
challenges lie ahead for the Operations Research community. Whereas the chemical 
shipper has about 4,000 lanes with 10 load locations to auction off, CPLEX can solve 
instances with no more than only 270 lanes with ten load locations to optimally. 

However, our simple randomized heuristic gives surprisingly good results, with 
the average solution quality no more than 0.8% away from optimal. Better per- 
formance may be expected from a more sophisticated heuristic, for example ran- 
domized local search. Hence, for practical purposes, the design of such randomized 
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heuristics seems to be a very promising avenue for further research. 

In our case, the shipper disallowed XOR bids. However, in practice, the car- 
riers value a bundle of lanes not always as the sum of the values of the individual 
lanes. A possible research direction is the analysis of a model with XOR-bids. The 
winner determination problem with XOR bids has been extensively analyzed (Sand- 
holm (2002) and Gtinluk et al. (2002)), without the much complicating constraints 
specific for the transportation sector, however. 
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LP * objective value of the optimal linear programming solution ~ 

ALH average gap between the best heuristic solution and LP* 

MLH maximal gap between the best heuristic solution and LP* 

Affl average gap between the best heuristic solution and the optimal integer solution 

j rnaximal * a P between the b ^ heuristic solution and the optimal integer solution 

I instances out of 25 with integer optimal LP solution 

#HO instances out of 25 with heuristic solution equal to optimal solution 

AT average total computation time in seconds 

MT maximal total computation time in seconds 

AL average computation time for solving LP relaxation 

ML rnaximal computation time for solving LP relaxation 

AN average number of nodes in the branch-and-bound tree 

MN maxima] number of nodes in the branch-and-bound tree 



Table 1: Notation used in reporting the computational results 
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Table 2: Computational results for the case MC=3, MC=4 and MC=5 
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Table 3: Computational results for the case MC=3, MC=4 and MC=. 
initial solutions 
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